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Abstract. Cosmology is an attracting subject for students but usually difficult to 
deal with if general relativity is not known. In this article, we first recall the Newtonian 
derivation of the Friedmann equations which govern the dynamics of our universe and 
discuss the validity of such a derivation. We then study the equations of evolution of 
the universe in terms of a dynamical system. This sums up the different behaviors of 
our universe and enables to address some cosmological problems. 

PACS numbers: ??? 
1. Introduction 

In this article, we want to present a pedagogical approach to the equations governing the 
evolution of the universe, namely the Friedmann equations. Indeed, the derivation of 
this equations is intrinsically relativistic. Although in Newtonian theory, the universe 
must be static, Milne Q and McCrea and Milne || showed that, surprisingly, the 
Friedmann equations can be derived from the simpler Newtonian theory. In section 
we recall their derivation (§ |2.1| ) for a universe filled with pressureless matter and then 
discuss the introduction of a cosmological constant (§ |2]^). Indeed, it is puzzling that 
the Newtonian theory and general relativity give the same results; we briefly discuss 
this issue in § EO[ . 

Once we have interpreted the Friedmann equations, we study them as a dynamical 
system. The first authors to consider such an approach were Stabell and Refsdal || 
who investigated the Friedmann-Lemaitre model with a pressureless fluid. This was 
then generalised to a fluid with any equation of state ||, |fl. Then, this technique 
was intensively used to study the isotropisation of homogeneous models (see e.g. |J 
and references therein). For a general description of the use of dynamical systems in 
cosmology, we refer to the book by Wainwright and Ellis [[!]] where most of the techniques 
are detailed. Our purpose here, is to present such an analysis for a fluid with any 
equation of state and including a cosmological constant while staying as pedagogical 
as possible. In section § [3], we rewrite the Friedmann equations under a form easier 
to handle with and we extract the dynamical system to study. We then determine the 
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fixed points of this system and discuss their stability. We illustrate this analytic study 
by a numerical integration of this set of equations (§ |]) and finish by a discussion about 
the initial conditions explaining the current observed state of our universe (§ |^). 

2. A Newtonian derivation of the Friedmann equation 

We follow the approach by Milne [I] and McCrea and Milne |2j and the reader is referred 
to H for further details. 

2.1. General derivation 

We consider a sphere of radius R filled with a pressureless fluid (P = 0) of uniform 
(mass) density p free-falling under its own gravitational field in an otherwise empty 
Euclidean space. We decompose the coordinate x of any particle of the fluid as 

x = a(i)r (1) 

where r is a constant vector referred to as the comoving coordinate, t is the time 
coordinate and a the scale factor. We choose a to have the dimension of a length 
and r to be dimensionless. It implies that the sphere undergoes a self similar expansion 
or contraction and that no particle can cross another one. Indeed the edge of the sphere 
is also moving as 

R(t) = a(t)R . (2) 

Assume that while sitting on a particle labelled % we are observing a particle labelled j; 
we see it drift with the relative velocity 

Vij = a (tj - Ti) = Hxij (3) 

where a dot refers to a time derivative, H = a/a and x^- = (r^ — r»). As a consequence, 
any particle i sees any other particle j with a radial velocity proportional to its distance 
and the expansion is isotropic with respect to any point of the sphere, whatever the 
function a(t). But, note that this does not imply that all particles are equivalent (as 
will be discussed later). 

To determine the equation of motion of any particle of this expanding sphere, we 
first write the equation of matter conservation stating that the mass within any comoving 
volume is constant (i.e. px 3 oc r 3 ) implying that 

pit) oc a~ 3 (t), (4) 

which can also be written under the form 

p + 3Hp = 0. (5) 

Note that Eq. (|) can also be deduced from the more general conservation equation 
dtp + Vaj = with j = pv, v = ifx and V x x = 3. 
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To determine the equation of evolution of the scale factor a, we first compute the 
gravitational potential energy Eq of a particle of masse m by applying the Gauss law 

Ea = - GM(< *> m (6) 

X 

where G is the Newton constant and M(< x) the mass within the sphere of radius x 
given by 

4:71 

M(< x) = —px 3 . (7) 

We then need to evaluate its kinetic energy Ek which takes the simple form 

E K = X -mx 2 . (8) 

The conservation of the total energy E = Eg + Ek implies, after the use of the 
decomposition (|J) and a simplification by r, that 

d\ 2 8ttG Kc 2 

P-—T (9) 



s.a) 3 o? 

where if is a dimensionless constant (which can depend on r) given by K = 
-2£/(mc 2 r 2 )|. 

2.2. Introducing a cosmological constant 

In the former derivation, the gravitational potential on any particle inside the sphere is 
proportional to the distance x 2 . Any other force deriving from a potential proportional 
to x 2 will mimic a gravitational effect. A force deriving from the potential energy E\ 
defined by 

Ac 2 

Ek = -m^x 2 (10) 
o 

where A is a constant was introduced by Einstein in 1917. As in the previous section, 
writing that the total energy E = E K + E G + E\ is constant leads to the equation of 
motion 

d\ 2 8ttG Kc 2 Ac 2 



-P- — + — . (11) 
,a/ 3 cr 3 

From (|To"D , we deduce that A has the dimension of an inverse squared length. The total 

force on a particle is 

F = ™ -— P + y Jx (12) 

from which it can be concluded that (i) it opposes gravity if A is positive and that (ii) 
it can be tuned so that F = leading to a = and p =constant if 

A = ^P- (13) 

% This scaling of K with r is imposed by the requirement that the expansion is self-similar (Eq. |l|) and 
that no shell of labeled r can cross a shell of label r' > r. 
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Table 1. Units of the quantities introduced in the article. M, L and T stand 
respectively for mass, length and time units. 

arv p P E H A K 

L - L.T- 1 ALL- 3 M.L-\T- 2 M.L 2 .T~ 2 T-l L~ 2 - 



This enables to recover a static autogravitating sphere hence leading to a model for a 
static universe. The force deriving from E\ is analogous to the one exerted by a spring 
of negative constant. 

To finish, we recall on table [lj the dimension of all the quantities used in the former 
sections, mainly to compare with standard textbooks in which the choice c = 1 is usually 
made. 

2.3. Discussion 

From this Newtonian approach, the equation of evolution of the universe identified 
with this gravitating sphere are thus given by equation (|) and flTTp. These are two 
differential equations for the two variables a(t) and p(t) which can be solved once the 
two parameters K and A have been chosen. 

In the context of general relativity, one can deduce the law of evolution for the scale 
factor of the universe a which is given by the Friedmann equations 

tt2 k Kc 2 Ac 2 
3 a z 3 

a k . P. Ac 2 , . 

- =-6<" + 3 ?» + - (15) 
with k = 8ttG and the conservation equation 

p + 3H(p + -) = 0. (16) 

Eq. (0) reduces to (0) and, Eq. ([TBD to (|) when P = 0. The equation (|15D is 
redundant and can be deduced from the two others. Note that now Eq. (|16D is also a 
conservation equation but with the mass flux j = (p + P/c 2 )v. This can be interpreted 
by remembering that the first law of thermodynamics for an adiabatic system takes 
the form 

E + PV = (17) 

where E = pVc 2 is the energy contained in the physical volume V (scaling as a 3 ). 

The first thing to stress is that equations (|5]) and (IT) do not depend on the radius 



Ro of the sphere. It thus seems that we can let it go to infinity without changing the 
conclusions and hence explaining why we recover the Friedmann equations. This was 
the viewpoint adopted by Milne |TJ and McCrea and Milne 0. This approach leads to 
some problems. First, it has to be checked that the Gauss theorem still applies after 
taking the limit toward infinity (i.e. one has to check that the integrals and the limit 
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commute). This imposes that p decreases fast enough with r and thus that there is 
indeed a center. Equivalently, as pointed out by Layzer flOfl , the force on any particle 
of an infinite homogeneous distribution is undetermined (the integral over the angles 
is zero while the integral over the radial coordinate is infinite). The convergence of 
the force requires either the mass distribution to be finite (in which case it can be 
homogeneous) or to be inhomogeneous if it is infinite. The issue of the finiteness of the 
universe has been widely discussed and a clear presentation of the evolution of ideas in 
that respect are presented in fl6|]. Second, for distances of cosmological interests, i.e. of 
some hundred of Megaparsec, the recession speed of the particles of the sphere are of 
order of some fraction of the speed of light. One will thus require a (special) relativistic 
treatment of the expanding sphere. Third, the gravitational potential grows with the 
square of the radius of the sphere but it can not become too large otherwise, due to the 
virial theorem, the velocities would exceed the speed of light. 

It was then proposed |L2| that such an expanding sphere may describe a region 
the size of which is small compared with the size of the observable universe (i.e. of the 
Hubble size). Since all regions of a uniform and isotropic universe expand the same way, 
the study of a small region gives information about the whole universe (but this does 
not solve the problem of the computation of the gravitational force). 

The center seems to be a privileged points since it is the only point to be at rest 
with respect to the absolute frame. But, one can show that the spacetime background of 
Newtonian mechanics is invariant under a larger group than the traditionally described 
Galilean group. As shown by Milne pj, McCrea and Milne |2| and Bonnor (see also 
Carter and Gaffet |TJ| for a modern description) it includes the set of all time-dependent 
space translations 

x l ^x i + z i (t) 

where z l (t) are arbitrarily differentiable functions depending only on the time coordinate 
t. This group of transformation is intermediate between the Galilean group and the 
group of all diffeomorphisms under which the Einstein theory in invariant. Thanks to 
this invariance group, each point can be chosen as a center around which there is local 
isotropy and homogeneity but the isotropy is broken by the existence of the boundary 
of the sphere (i.e. all observer can believe living at the center as long as he/she does 
not observe the boundary of the expanding sphere). 

There are also conceptual differences between the Newtonian cosmology and the 
relativist cosmology. In the former we have a sphere of particle moving in a static and 
absolute Euclidean space and the time of evolution of the sphere is disconnected from 
the absolute time t. For instance in a recollapsing sphere, the time will go on flowing 
even after the crunch of the sphere. In general relativity, space is expanding and the 
particles are comoving. We thus identify an expanding sphere in a fixed background 
and an expanding spacetime with fixed particles. As long as we are dealing with a 
pressureless fluid, this is possible since there is no pressure gradient and each point of 
the sphere can be identify with one point of space (in fact, with an absolute time we are 
working in a synchronous reference frame and we want it to be also comoving, which is 
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Table 2. Comparison of the nature of the Newtonian trajectory and of the structure 
of space according to the value of the constant K in Eq. (fnl) . 



E 


> 





< 


Trajectory 


hyperbolic 
unbounded 


parabolic 
unbounded 


elliptic 
bounded 


K 


< 





> 


Spatial section 


infinite 


infinite 


finite 



possible only if P = flI5f). Moreover, the pressure term in the Friedmann equations 
cannot trivially be recovered from the Newtonian argument. As shown, one gets the 
correct Friedmann equations if one starts from the conservation law including pressure 
(and derived from the first law of thermodynamics) and the conservation of energy But 
if one were starting from the Newton law relating force (|12D and acceleration (max), the 
term containing the pressure in fll5|) would not have been recovered; one should have 
added an extra pressure contribution Fp = —AnGmPit/c 2 which can not be guessed. 
This is a consequence that in general relativity any type of energy has a gravitational 
effect. In a way it is a "miracle" that the equation (14) does not depend on P, which 
makes it possible to derive from the Newtonian conservation of energy. Beside it has also 
to be stressed that the Newtonian derivation of the Friedmann equations by Milne came 
after Friedmann and Lemaitre demonstrated the validity of the Friedmann equations 
for an unbounded homogeneous distribution of matter (using general relativity). It has 
to be pointed out that these Newtonian models can not explain all the observational 
relations since, contrary to general relativity, they do not incorporate a theory of light 
propagation. As outlined by Lazer [TIJ one can sometime legitimately treat a part 
of the (dust) expanding universe as an isolated system in which case the Newtonian 
treatment is correct, which makes McCrea JTT) conclude that this is an indication that 
Einstein's law of gravity must admit the same interpretation as that of Newton's in the 
case of a spherically symmetric mass distribution. Note that the structural similarity of 
Einstein and Newton gravity were put forward by Cartan [17| who showed that these 
two theories are much closer that one naively thought and, in that framework (which 
goes far beyond our purpose) one can work out a correct derivation of the Friedmann 
equations (see e.g. @)- 

The most important outcome of the Newtonian derivation of the Friedmann 
equations is that it allows to interpret equation (|HD in terms of the conservation of 
energy; the term in H 2 represents the kinetic energy, the term in np/3 the gravitational 
potential energy, the term in A/3 the energy associated with the cosmo logical constant 
and the term in K the total energy of the system. The properties of the spatial sections 
(i.e. of the three dimensional spaces of constant time) are related to the sign of K and 
can be compared with the property of the trajectories of the point of the sphere which 
are related to the sign of the total energy E; we sum up all these properties on table 0. 
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3. The Friedmann equations as a dynamical system 

The Friedmann equations (|T^-|1^) and the conservation equation ([16|) form a set of two 
independent equations for three variables (a, P and p). The usual approach is to solve 
this system by specifying the matter content of the universe mainly by assuming an 
equation of state of the form 

P = (7 - l)pc 2 (18) 

where 7 may depend on p and thus on time. For a pressureless fluid (modelling 
for instance a fluid of galaxies) 7 = 1 and for a fluid of radiation (such as photon, 
neutrino,...) 7 = 4/3. We assume that 7 7^ since such a type of matter is described 
by the cosmological constant and singled out from "ordinary" matter and that 7 7^ 2/3 
since such a type of matter mimics the curvature term and is thus incorporated with it. 

One can then first integrate ( |i~6l) rewritten as dp/ (7/)) = — 3da/a to get the function 
p(a) which, in the case where 7 is constant, yields 

p(a) = Co -37 (19) 

where C is a positive constant of integration, and then insert the solution for p(a) in 
Eq. fll4]) to get a closed equation for the scale factor a (see e.g. |0| for such an approach 
and |20| for an alternative and pedagogical derivation). 

In this section, we want to present another approach in which the Friedmann 
equations are considered as a dynamical system and to determine its phase space. 

3.1. Derivation of the system 

The first step is to rewrite the set of dynamical equations with the three new variables 
Q, Q\ and Q# defined as 

" < 20 » 

Ac 2 

Q* = m . (2D 
Kc 2 

n *=-Hnk- ( 22 ) 

They respectively represent the relative amount of energy density present in the matter 
distribution, cosmological constant and curvature. Q has to be positive and there is no 
constraint on the sign of both and Qk- With these definitions, it is straightforward 
to deduce from (|TJ) that 

n + n A + n K = 1. (23) 

Using that H = d/a — H 2 , expressing d/a from Eq. (|T5| ) and H 2 from Eq. (|14"D , we 
deduce that 

|s = -(i + g) (24) 



A Dynamical Study of the Friedmann Equations 



8 



where the deceleration parameter q is defined by 



q 



3^ 



37 



si, 



(25) 



2 v 7 2 

It is useful to rewrite the full set of equations by introducing the new dimensionless time 
variable rj = ln(a/ao), ao being for instance the value of a today. The derivative of any 
quantity X with respect to r], X', is then related to its derivative with respect to t by 
X' — X/H. The equation of evolution of the Hubble parameter (]24 ) takes the form 

H' = -(l + q)H. (26) 

Now, differentiating Q, Q A and Q K with respect to 77, using Eq. ( f26|) to express H', 
a' = a and Eq. ([U|) to express p' = — 37p, we obtain the system 

fi' = (2g + 2 - 37)^ 



(27) 
(28) 
(29) 



Indeed, it is useless to study the full system (|2^-^D|) (i) since H does not enter the 
set of equations (|27^p9|) and is solely determined by Eq. ( pB"D once this system has been 
solved and (ii) since fl can be deduced algebraically from (E3p. As a consequence, we 
retain the closed system 



n' A = 2(l + q)Sl A 
Q' K = 2qQ K 

and it is trivial to check that Q f + + Q' K = as expected form 



(30) 



f Q' A = 2{l + q)Q A 

\n> K = 2 g n K 

with q being a function of fl A and Qk only and defined in (f25|). 

The system ( PD|) is autonomous |^T|, which implies that there is a unique integral 
curve passing through a given point, except where the tangent vector is not defined 
(fixed points). Note that at every point on the curve the system ( |3~0"1) assigns a unique 
tangent vector to the curve at that point. It immediately follows that two trajectories 



cannot cross; otherwise the tangent vector at the crossing point would not be unique 21 



3.2. Determination of the fixed points 

To study the system ( |30D as a dynamical system, we first need to determine the set of 
fixed points, i.e. the set of solutions such that Q' A = and Q' K = 0. These solutions 
represent equilibrium positions which indeed can be either stable or unstable. The fixed 
points are thus solutions of 

(1 + q )n A = 0, qVl K = 0. (31) 

We obtain the three solutions 



(0^,00 e {(0,0), (0,1), (1,0)}. (32) 

Each of these solutions represent a universe with different physical characteristics: 
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(i) (fl K ,Q A ) = (0,0): the Einstein de Sitter space (EdS). 

It is a universe with flat spatial sections, i.e. the three dimensional hypersurfaces of 
constant time are Euclidean and it has no cosmological constant. We deduce from 
([2D and (f2|) that 

Q = 1, q = ~ 7 - 1 (33) 
and integrating Eq. (|14D gives 

«(*) = (^X*) (34) 
for the solution vanishing at £ = 0. 

(ii) (Qk,Qa) = (0, 1): the de Sitter space (dS). 

It is an empty space filled with a positive cosmological constant and with flat spatial 
sections. We deduce from (E3T) and (1251) that 



Q = o, g = -1 (35) 
and integrating Eq. (|l4l) gives 

a(i) = aoev 7 ?*. (36) 
This universe is accelerating in an eternal exponential expansion. 

(iii) (Qk,&>a) = (1)0): the Milne universe (M). 

It is an empty space with no cosmological constant and with hyperbolic spatial 
section (K < 0). We deduce from @ and (|2|) that 

= 0, q = (37) 

and integrating Eq. (|14|) gives 

a(t) = dot (38) 

for the solution vanishing at t = and in units where X = tig. 

It is also interesting to study the properties of the three following invariant lines 
which separate the phase space in disconnected regions: 

(i) VLk = : The system ([3(]) reduces to the equation of evolution for fl\ 

n' A = 3 7 (1 - Q A )tt A . (39) 

Thus, if initially Qk — 0, we stay on this line during the whole evolution and 
converge toward either £7a = 1 (i-e. the fixed point dS) or toward = (i.e. the 
fixed point EdS). It also follows that no integral flow lines of the system ([Jl]) can 
cross the line Qk = 0. It separates the universes with Qk > which are compact 
(i.e. having a finite spatial extension) and the universes with VLk < which are 
infinite (if one assumes trivial topology P3|). Crossing the line Qk — would thus 
imply a change of topology. Note that if 7 = 0, the fluid behaves like a cosmological 
constant and thus Q' A = since Qk remains zero. 
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(ii) Q A = : The system fl3~Q|) reduces to the equation of evolution for Q K 

Q' K = (3 1 -2)(1~Q K )Q K . (40) 

As in the previous case, we stay on this line during the whole evolution and converge 
toward either — 1 (i-e. the fixed point M) or toward ^a = (i.e. the fixed 
point EdS). It also follows that no integral flow lines of the system (^) can cross 
the line Q A = 0. Note that if 7 = 2/3, the fluid behaves like a curvature term and 
thus Q' K = since Q A remains zero. 

(iii) Q = : It is a boundary of the phase space since Q is non negative. We now have 
q = —Q A and the system (|3~0"D reduces to 

n' A = 2(1 - n A )n A . (4i) 

The universe converges either toward (dS) or (M). 



3.3. Stability analysis 

The second step is to determine whether these fixed points are stable (i.e. attractors: 
A), unstable (i.e. repulsor: R) or saddle (S) points (i.e. attractor in a direction and 
repulsor in another). This property can be obtained by studying at the evolution of a 
small deviation from the equilibrium configuration. We thus decompose fix and Q A as 

VL K = VL k + u) K , (42) 
Q A =U A + u A , (43) 

where (fl A , Qk) represents the coordinates of one of the fixed points determined in the 
previous section and where (u a ,uk) is a small deviation around this point. 
Writing the system of evolution (|30"D as 

( n K \ = ( F K (n A ,n K ) \ 

\n A J { F A (n A ,n K ) ) ' 1 ] 

where Fr and Fa are two functions determined from (|30D , it can be expanded to linear 
order around (Qk,^a) (for which Fx and Fa vanish) to give the equation of evolution 
of (u A ,u K ) 

dFjL 

" ■„ B) l£J-"*™u:j- (45) 

The stability of a given fixed point depends on the sign of the two eigenvalues (Ai^) of the 
matrix P(q a Q K y If both eigenvalues are positive (resp. negative) then the fixed point is a 
repulsor (resp. an attractor) since (ujk,w a ) will respectively goes to infinity (resp. zero). 
In the case where the two eigenvalues have different signs, the fixed point is an attractor 
along the direction of the eigenvector associated with the negative eigenvalue and a 
repulsor along the direction of the eigenvector associated with the positive eigenvalue. 
We also introduce u Al 2 the eigenvectors associated to the two eigenvalues which give 
the (eigen)-directions of attraction or repulsion. 





A Dynamical Study of the Friedmann Equations 



11 



We have to perform this stability analysis for each of the three fixed points 
(reminding that 7 7^ 0, 2/3): 

(i) EdS fixed point: In that case, the matrix V is given by 

M 3 v 2 3°> (46) 

the eigenvalues of which are trivially given by 37 — 2 and 27. We thus conclude that 
if 7 G] — 00, 0[ then EdS is an attractor, that it is a saddle point when 7 g]0, 2/3[ and 
a repulsor when 7 G]2/3, +00 [. The matrix Ve&s being diagonal the two eigenvectors 
are trivially given by 

u ( 3 7) = (M), U ( 3 7 _ 2 ) = (1,0) (47) 
corresponding respectively to two invariant boundaries Q\ = and Qk = 0. 

(ii) dS fixed point: The matrix V is now given by 

^=( 2 : 2 37 X)- m 

The eigenvalues of Vds are — 2 and —37. It follows that the fixed point dS is never 
a repulsor. If 7 g] — 00, 0[ then dS is a saddle point and, when 7 g]0, +00 [, it is an 
attractor. The two eigenvectors are now given by 

U ( _3 7 ) = (0,1), U(_ 2) = (l,-1) (49) 

corresponding respectively to the two boundaries Q = and = 0. 

(iii) M fixed point: The matrix V is now given by 



Pm= „ ' T ■ (60) 




The eigenvalues of Vm are 2 and 2 — 37. It follows that M is never an attractor 
since one of his eigenvalues is always positive. If 7 g] — 00, 2/3[ then M is a repulsor 
point and, when 7 G]2/3, +00 [, it is a saddle point. The two eigenvectors are now 
given by 

Up-atf = (1,0), U( 2 ) = (1,-1) (51) 

corresponding respectively to the two boundaries Qk — and O = 0. 

Before we sum up all theses results, let us concentrate about the cases where 7 = 
or 7 = 2/3 in which the matter behaves respectively either as a cosmo logical constant 
or as a curvature term. As a consequence Q can be absorbed in a redefinition of either 
Q\ or Qk and we can set Q = from which it follows that (|23| ) implies + = 1- 
In both cases, we deduce from ([250 that q = Qk — 1 = — so that 

tt' K = 2(tt K - l)tt K , tt' A = 2(1 - tt A )tt A (52) 

which are not independent equations due to the constraint + = 1- Thus, for 
7 = or 7 = 2/3, the two fixed points are either (M) or (dS) which are respectively a 
repulsor and an attractor. 
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Table 3. Stability properties of the three fixed points (EdS, dS and M) as a function 
of the polytropic index 7. (A: attractor, R: repulsor and S: saddle point) 
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Figure 1. The fixed points and their stability depending of the value of the index 7: 
(a) 7 < 0, (b) < 7 < 2/3 and (c) 7 > 2/3. 



As a conclusion of this study, we sum up the properties of the three spacetimes 
as a function of the polytropic index of the cosmic fluid in table ^| and in figure |l], 
we depict the fixed points, their directions of stability and instability as well as the 
invariant boundary in the plane &k)- Indeed, the attractor solution can be guessed 
directly from Eq. (|14D and the behavior ( ]19|) of the density with the scale factor since 
if 7 < the matter energy density scales as a -37 and comes to dominate over the 
cosmological constant (scaling as a ) and the curvature (scaling as a -2 ). On the other 
hand the cosmological constant always finishes by dominating if 7 > 0. The curvature 
can never dominates in the long run since it will be caught up by either the matter or 
the cosmological constant. 



4. Numerical examples 

The full phase space picture can be obtained only through a numerical integration of 
the system ( j30|) by using an implicit fourth order Runge-Kutta method . 



Ordinary matter such as a pressureless fluid or a radiation fluid has 7 > 1 and we 
first consider this case on figure ^| where we depict the phase space both in the (fix, ^a) 
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Figure 2. Phase space of the system ( |30| ) in the plane (Qk,^a) [left] and (f2, 
[right] . We have represented the three fixed points and the lines Qk — 1 and 17a = 1 
and we have considered the value 7 = 1 (i.e. pressureless fluid). 
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Figure 3. Phase space of the system ( p0| ) in the plane (Qk,^a) [left] an d ^a) 
[right] for 7 = 1/3. 



where the analytic study of the fixed points was performed but also in the plane (Q, 
for complementarity. On figure we consider the case where < 7 < 2/3 which can 
corresponds to a scalar field slowly rolling down its potential or a tangle of domain 
strings (for which 7 = 1/3) and we finish by the more theoretical case where 7 < on 
figure |] for which we know no simple physical example (see however 
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5. Discussion and conclusions 



To discuss the naturalness of the initial conditions leading to our observed universe, we 
have to add the actual observational measures in the plane (Q, Oa) and trace them back 
to estimate the domain in which our universe has started. This required (i) to know 
what are the constraints on the cosmological constant and the curvature of the universe 
and (ii) determine the age of the universe, i.e. the time during which we must integrate 
back. 

It is not the purpose of this article to detail the observational methods used in 



cosmology and a description can be found en e.g. ||19|| ; we now just sum up what is 
thought to be the current status of these observations. The current observational data 
such as the cosmic microwave background measurements El, the Type la supernovae 



data large scale velocity fields |26|| , gravitational lensing |27J and the measure of 



the mass to light ratio |28| tend to show that 

tt ~ 0.3, tt Ao ~ 0.7. (53) 



We refer the reader to the review by Bahcall et al. [2§] for a combined study of these 
data and a description of all the observation methods. Let us just keep in mind that we 
are close to the line fix — and let us consider the safe area of parameter such that 

V : {n E [0.1,0.5], Ov G [0.5,0.9]} (54) 

and let us determinate the initial conditions allowed by these observations. 

For that purpose, we need to integrate the system ( |30"D back in time during a time 
equal to the age of the universe. Today, the matter content of the universe is dominated 
by a pressureless fluid, the energy density of which is obtained once f2 has been chosen 
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and is 

p ma t = (-) ~ 3 = 1-80 x 10 - 29 fi (-) ~ 3 g.cm- 3 (55) 

K \a / \a / 

where H = 100 h km/s/Mpc is the Hubble constant today. The energy density of 

the radiation is obtained by computing the energy contained in the cosmic microwave 

background which is the dominant contribution to the radiation in the universe. Since it 

is a black body with temperature 6 = 2.726 K, we deduce, from the Stephan-Boltzmann 

law, that 

-34 ( a \ 4 „ _-3 



p rad = 4.47(1 + /„) x 10" 34 - g.cm" 3 (56) 

\a / 

where f v = 0.68 is a factor to take into account the contribution of three families of 
neutrinos [TP]. The radiation was thus dominating over the matter for scale factors 
smaller than a eq at which p mat = p rad and thus given by 



a eq 4.5 x 10 



-5 



ao \l n z 

We can integrate back until the Planck era for which ctpi/ao ~ 10~ 30 and can thus 
approximate 7 by 

fl a eq <a<a 
[ 4/3 api < a < a eq 

which is a good approximation for 7 = 1 + 1/3(1 + a/a eq ). In figure ||, we depict the 
domain T>q of current observational values and its inverse image by the system (|30]) at 
the beginning of the matter era and at the end Planck era. 

To illustrate this fine tuning problem analytically, let us just consider the simplest 
case where fl\ = for which the evolution of fl is simply given by 

Q' = (37 - 2)Q(Q - 1) (59) 



the solution of which is 

Q V ao J 



n = o.. , ,37-2 (6°) 



and thus, 



S2q a 



< a < a 



Q = { T . . .. (61) 



From the observational data we get that Qk ~ O(10 x ) and thus Qo ~ from 
which we deduce that 

^U cq ~ 0(1O" 4 ), n*|«=« PI ~ O(10^). (62) 

This illustrate that a almost flat universe requires a fine tuning of the curvature, which 
is a consequence that = is a repulsor in both the radiation and the matter era. 
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0,2 0,4 0,6 0,8 1 0,2 0,4 0,6 0,8 1 

12 12 

A A 

Figure 5. The domain T>o of current observational values for our universe and its 
inverse image by the system (|30| ) when one is assuming that the universe has been 
dominated by a pressureless fluid during its whole evolution. 

A solution to solve this fine tuning problem would be to add a phase prior to the 
radiation era in which 7 < so that Qk — becomes an attractor and to tune the 
duration of this era such as to have the correct initial conditions. Then, during the 
standard evolution EdS becomes repulsor and we evolve toward dS staying close to the 
line Qk — hence explaining current observations. Inflation is a realisation of such 
a scenario. What inflation really does is to change the stability property of the fixed 
points and invariant boundaries of the system (|30|). Hence during this period where a 
fluid with negative pressure is dominating we are attracted close to the line fix = and 
the closer the longer this phase lasts. We then switch to a phase with normal matter and 
start to drift away due to repulsive property of EdS. Nevertheless, inflation does more 
than just explaining where we stand in this phase space, it also gives an explanation 
for the observed structures (galaxies, clusters...) of our universe, but this is beyond the 
scope of this article. 
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